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ABSTRACT 

This  paper  considers  scattering  screens  that  have  arbitrary  spatial  variations  of  scattering  strength 
transverse  to  the  line  of  sight,  including  screens  that  are  spatially  well  confined,  such  as  disks  and  fila¬ 
ments.  We  calculate  the  scattered  image  of  a  point  source  and  the  observed  pulse  shape  of  a  scattered 
impulse.  The  consequences  of  screen  confinement  include  (1)  source  image  shapes  that  are  determined  by 
the  physical  extent  of  the  screen  rather  than  by  the  shapes  of  much  smaller  diffracting  microirregularities 
(these  include  image  elongations  and  orientations  that  are  frequency  dependent);  (2)  variation  with  fre¬ 
quency  of  angular  broadening  that  is  much  weaker  than  the  trademark  v~2  scaling  law  (for  a  cold, 
unmagnetized  plasma),  including  frequency-independent  cases;  and  (3)  similar  departure  of  the  pulse¬ 
broadening  time  from  the  usually  expected  v-4  scaling  law.  We  briefly  discuss  applications  that  include 
scattering  of  pulses  from  the  Crab  pulsar  by  filaments  in  the  Crab  Nebula;  image  asymmetries  from 
Galactic  scattering  of  the  sources  Cyg  X-3,  Sgr  A*,  and  NGC  6334B;  and  scattering  of  background 
active  galactic  nuclei  by  intervening  galaxies.  We  also  address  the  consequences  for  inferences  about  the 
shape  of  the  wavenumber  spectrum  of  electron  density  irregularities,  which  depend  on  scaling  laws  for 
the  image  size  and  the  pulse  broadening.  Future  low-frequency  (<100  MHz)  array  observations  will  also 
be  strongly  affected  by  the  Galactic  structure  of  scattering  material.  Our  formalism  is  derived  in  the 
context  of  radio  scattering  by  plasma  density  fluctuations.  It  is  also  applicable  to  optical,  UV,  and  X-ray 
scattering  by  grains  in  the  interstellar  medium. 

Subject  headings :  ISM :  structure  —  pulsars :  general  —  scattering 


1.  INTRODUCTION 

Images  of  scattered  radio  sources  and  distorted  pulses  from  pulsars  provide  some  of  the  most  used  observables  for  probing 
microstructure  in  the  electron  density  of  interstellar  gas.  Over  the  last  decade,  interstellar  scattering  measurements  have 
revealed  asymmetries  in  the  scattered  images  of  radio  sources.  These  are  most  often  interpreted  in  terms  of  underlying 
anisotropy  of  the  very  small  irregularities  that  diffract  the  radiation.  That  anisotropy,  in  turn,  most  likely  reflects  the 
orientation  of  magnetic  fields  in  the  H  n  gas  that  contains  the  microstructure.  Angular  broadening  of  compact  sources  and 
pulse  distortions  due  to  multipath  propagation  are  used  to  probe  the  amplitude  of  scattering  and  also,  through  the  frequency 
scaling,  to  constrain  the  shape  of  the  wavenumber  spectrum  for  the  microstructure.  Inversion  of  scattering  observables  into 
information  about  the  micro  structure  almost  invariably  relies  on  the  assumption  that  the  scattering  strength  is  uniform  in 
directions  transverse  to  the  line  of  sight. 

We  reconsider  the  assumptions  used  to  analyze  angular  and  temporal  broadening,  in  particular  the  assumption  of 
uniformity  of  the  scattering  medium  transverse  to  the  line  of  sight.  One  reason  is  that  the  interstellar  medium  (ISM)  shows 
structures  on  a  wide  variety  of  scales  and  so  it  is  reasonable  to  expect  manifestations  of  nonuniformities,  at  least  in  some 
directions.  Second,  the  physics  that  underlies  asymmetric  images  is  quite  different  if  the  asymmetry  occurs  on  scales  much 
larger  than  diffractive  scales,  as  they  would  if  the  asymmetry  is  caused  by  the  large-scale  distribution  of  diffracting  irregu¬ 
larities.  Third,  observations  of  the  Crab  pulsar  show  anomalous  scalings  of  pulse  broadening  with  frequency.  These  are 
interpreted  by  some  as  indicating  that  scattering  occurs  within  the  pulsar  magnetosphere  rather  than  in  a  cold  plasma  (J. 
Eilek  1997,  private  communication;  Hankins  &  Moffett  1998;  Lyutikov  &  Parikh  2000).  As  we  show,  anomalous  scalings 
occur  quite  naturally  from  cold  plasma  extrinsic  to  the  pulsar  if  the  scattering  region  is  bounded  in  the  transverse  direction. 
Additionally,  future  observations  at  low  radio  frequencies  of  a  variety  of  sources — including  high-redshift  sources — are 
expected  to  reveal  further  anomalous  scattering  that  most  likely  will  be  the  result  of  confined  scattering  structures.  Finally,  the 
scaling  with  frequency  of  angular  and  pulse  broadening  is  often  used  to  constrain  the  shape  of  the  wavenumber  spectrum  of 
scattering  irregularities  (e.g.,  Cordes,  Weisberg,  &  Boriakoff  1985 ;  Fey,  Spangler,  &  Cordes  1991).  Weakening  of  the  frequency 
dependence  by  confined  plasma  structures  would  be  interpreted  as  a  steeper  wavenumber  spectrum.  Thus,  it  is  important  to 
assess  the  role  of  confined  structures  in  the  observations  of  scattered  radio  sources. 

In  §  2  we  discuss  previous  treatments  of  angular  and  pulse  broadening.  In  §  3  we  derive  a  general  formalism  for  scattering 
that  takes  into  account  arbitrary  variations  of  scattering  strength  transverse  to  the  line  of  sight.  Examples  are  given  in  §  4. 
Applications  to  the  Crab  pulsar  and  other  Galactic  sources  are  given  in  §  5.  Future  observations  of  extragalactic  sources  and 
at  low  frequencies  are  considered  in  §  6.  In  §  7  we  discuss  possible  implications  for  the  interstellar  medium.  Finally,  in  §  8  we 
summarize  our  results. 
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2.  PAST  TREATMENTS  OF  ANGULAR  AND  PULSE  BROADENING 

It  is  well  known  that  the  shape  of  a  scattered  impulse,  viewed  through  a  thin  (along  the  line  of  sight),  infinitely  extended 
(transverse  to  the  line  of  sight)  screen  with  a  circularly  symmetric,  Gaussian  angle  distribution  is  a  one-sided  exponential 
function  (e.g.,  Rickett  1990  and  references  therein).  Thick  scattering  screens  produce  slower  rise  times,  while  screens  contain¬ 
ing  Kolmogorov  irregularities  (e.g.,  Rickett  1990;  Lambert  &  Rickett  1999)  produce  decays  that  are  slower  than  exponential. 
Results  along  these  lines  have  been  presented  by  Scheuer  (1968),  Williamson  (1972,  1973,  1975),  Lee  &  Jokipii  (1975),  and 
Isaacman  &  Rankin  (1977).  Williamson  (1975)  has  shown  that  the  pulse-broadening  function  from  multiple,  discrete  screens 
or  from  a  continuous  medium  consists  of  an  n-fold  convolution  of  (one  sided)  exponential  functions.  Williamson’s  result 
applies  to  the  case  where  the  phase  structure  function  is  square  law  in  form.  Media  with  Kolmorogov  wavenumber  spectra 
produce  different  shapes,  though  the  differences  are  small  compared  to  effects  that  we  consider  in  this  paper.  The  key,  implicit 
assumption  in  Williamson’s  analysis  (and  essentially  all  other  published  results  on  interstellar  pulse  broadening;  however,  see 
Lyne  &  Thorne  1975)  is  that  the  transverse  extent  of  the  scattering  screen  is  arbitrarily  large  and  that  the  strength  of  scattering 
is  uniform  across  the  screen.  If  these  assumptions  are  relaxed,  quite  different  results  emerge. 

We  show  that  several  “  anomalous  ”  phenomena  occur  when  scattering  structures  have  finite  transverse  extents,  including 

1.  Angular  broadening  from  scattering  or  refraction  that  scales  anomalously  with  v.  In  the  case  of  radio  scattering  in  cold 
plasmas,  anomalous  scaling  is  defined  as  a  significant  departure  from  a  v~ 2  scaling,  which  is  determined  by  the  microphysics 
of  the  plasma.1  Departures  will  always  consist  of  a  weaker  dependence  on  frequency. 

2.  Elongated  (or  otherwise  distorted)  images  of  point  sources  that  are  due  to  scattering  but  do  not  scale  as  v  " 2 . 

3.  Multiple  imaging  by  multiple,  discrete  screens  with  image  intensities  influenced  by  dilution  from  scattering. 

4.  Temporal  broadening  of  pulsar  pulses  which  shows  a  weaker  dependence  on  frequency  than  v~4,  in  accord  with  the 
scaling  of  angular  broadening.2 

5.  Replication  of  pulses  by  multiple  imaging  from  an  ensemble  of  screens. 

The  goal  of  this  paper  is  to  discuss  the  impact  of  confined  or  heavily  modulated  scattering  screens  on  some  of  the  basic 
observables,  primarily  pulse  broadening  and  angular  broadening.  Therefore,  we  do  not  explicitly  consider  refraction  from 
large-scale  features  in  the  ISM.  We  do  so  for  two  reasons.  First,  purely  diffractive  effects  are  rich  enough  in  variety  that  we 
need  to  isolate  the  discussion  to  those  effects.  Second,  our  results  can  also  be  applied  to  some  cases  where  refraction  is 
important  by  considering  a  “renormalized”  version  of  the  Kirchoff  diffraction  integral  (e.g.,  Cordes,  Pidwerbetsky,  & 
Lovelace  1986).  Renormalization  of  the  large-scale  gradient  and  phase  curvature  caused  by  refraction  at  a  given  screen 
location  can  be  cast  as  an  increased  or  decreased  image  intensity  and  also  as  a  change  in  the  ellipticity  of  the  angular 
distribution  of  scattered  radiation.  In  the  following  we  use  a  probability  density,/,,,  to  describe  the  scattering  angles  a  from  a 
screen.  The  renormalization  approach  allows  some  of  the  effects  of  refraction — image  shifts  and  shapes  and  intensity 
changes — to  be  absorbed  into  /„.  Our  approach  considers  only  geometrical  path  length  contributions  to  arrival  times  and 
excludes  dispersive  delays,  which  are  associated  with  the  screen  itself.  For  some  situations,  dispersive  delays  can  be  important. 
In  this  paper,  the  points  that  we  wish  to  make  concern  the  geometrical  phase  and  a  complete  discussion  that  includes 
dispersive  delays  would  distract  our  discussion  of  these  points.  Consequently,  we  defer  to  another  paper  a  complete  treatment 
that  includes  all  contributions  to  arrival  times. 

3.  PROBABILITY  DENSITIES  FOR  ANGLE  OF  ARRIVAL  AND  TIME  OF  ARRIVAL 

We  derive  the  image  of  a  scattered  point  source  and  the  scattered  pulse  shape  of  an  impulse  by  calculating,  respectively,  the 
probability  density  functions  (PDFs)  for  the  angle  of  arrival  (AO A)  and  the  time  of  arrival  (TO A).  In  the  following,  we 
calculate  the  effects  of  scattering  while  using  some  of  the  language  and  mathematics  of  ray  theory.  Williamson  (1975)  has 
shown  the  equivalence  of  wave  and  ray  optics  for  some  contexts,  as  have  Cordes  &  Rickett  (1998).  In  this  paper,  we  make  the 
simplifying  assumption  that  the  TO  A  is  related  geometrically  to  the  AO  A.  This  relationship  applies  when  only  the  geometri¬ 
cal  path  length  of  a  given  ray  path  contributes  significantly  to  the  TO  A.  In  general,  the  TO  A  includes  another  term  related  to 
the  integrated  refractive  index  (the  integrated  electron  density  in  the  case  of  a  cold  plasma).  We  ignore  the  nongeometrical 
term  because  there  are  astrophysical  contexts  where  its  contribution  is  negligible. 

Consider  a  series  of  diffracting  screens  at  distances  Dsj  from  a  source  that  is  at  distance  D  from  the  observer.  Letting  be 
the  (two  dimensional)  scattering  angle  from  the  jth  screen,  the  angular  deviation  0(s)  of  a  ray  path  and  its  transverse  offset  x(s) 
from  the  direct  ray  path  at  distance  s  from  the  source  are 


0(S)  =  0;  +  X  UJ  U(S  -  Dsj)  > 

j 

(1) 

x(s)  =  S0;  +  X  (s  -  Dsj)aj  U(s  -  Dsj) , 

(2) 

j 


where  U(x )  is  the  unit  step  function.  The  first  equation  relates  the  observed  ray  angle  ( 0 )  to  the  initial  ray  angle  (0t)  and  the 
scattering  angles  ( aj ).  The  relation  x(D)  =  0  stipulates  that  rays  must  reach  the  observer.  We  assume  that  all  angles  are  small 


1  Scattering  from  irregularities  with  a  Kolmogorov  wavenumber  spectrum  shows  v_u/5  scaling  under  the  circumstance  of  moderate  scattering  (see 
Cordes  &  Lazio  1991).  We  do  not  consider  this  scaling  anomalous.  When  a  confined  structure  contains  Kolmogorov  fluctuations,  the  scaling  of  angular  size 
may  be  considerably  shallower  than  v  " 1 1/5. 

2  A  Kolmogorov  spectrum  can  show  a  v  " 22/5  scaling  for  pulse  broadening.  As  with  angular  broadening,  we  do  not  consider  this  anomalous. 
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( 1 0t | ,  \cij\,  |  0 1  1),  though  it  is  not  difficult  to  extend  our  results  to  large  angles.  Note  that  cij  is  a  random  variable  described 

by  a  distribution  of  angles  that  is  determined  by  diffraction  (and,  as  mentioned  above,  can  also  include  refraction). 

Including  only  the  geometric  path  length  difference,  the  corresponding  time  delay  relative  to  the  direct  ray  path  is 

t  =  ^\  ds\0(s)\2  •  (3) 

The  overall  time  delay  also  includes  dispersive  components  which,  as  stated  above,  we  choose  to  ignore  because  our  main 
points  concern  the  effects  of  truncated  screens  on  the  arrival  times. 

To  calculate  the  probability  density  function  (PDF)  of  the  observed  angle  of  arrival,  6,  and  the  time  of  arrival,  t,  we  use 
Dirac  delta  functions  to  enforce  equation  (1)  and  x(D )  =  0  for  those  rays  that  reach  the  observer.  We  use  conditional 
probabilities  to  include  these  relations  and  to  integrate  over  the  PDFs  for  the  scattering  angles  in  each  screen  and  over  the 
PDF  of  0f.  The  result  is  simple  for  an  isotropic  source  or,  less  restrictively,  where  the  PDF  of  0t  is  constant  over  the  relevant 
range  of  initial  ray  angles,  6t,  as  we  assume.  Note  that  the  equation  x(D)  =  0  allows  us  to  eliminate  0t  as  an  independent 
variable. 

It  is  standard  to  assume  that  the  scattering  strength  is  invariant  across  a  scattering  screen.  Here  we  specify  a  more  general 
description.  Consider  each  screen  to  scatter  or  refract  rays  according  to  a  PDF  f  whose  width  varies  arbitrarily  across  the 
screen.  Accordingly,  we  write  each  screen’s  PDF  as /„  .[«,■;  x(D  •)],  where  x(D  )  is  equation  (2)  evaluated  at  the  location  of  each 
screen,  s  =  D 

Let  Q  be  an  observed  quantity  such  as  the  AOA  or  TOA  and  let  Q'(0t,  aj,  Dsj,  D )  be  its  value  given  0t,  aj}  and  Dsj  (and  it  is 
implicit  that  we  consider  all  j  when  there  are  multiple  screens).  The  PDF  of  Q  is 

,  =  EL  I  dajLlay,  xjD,^  j  -  Q'(0„  ajf  Dsj,  D)] 

2  n,  I  ^  x(Dsjn  ]'  dd;  /e/0i)<5[x(D)]  ’  U 

where  the  numerator  is  the  joint  PDF  of  Q  and  x(D)  =  0  (that  rays  reach  the  observer),  while  the  denominator  is  the  PDF  that 
rays  reach  the  observer.  Using  equation  (2),  we  transform  <5[(x(D)]  to  a  delta  function  involving  6t,  perform  the  integral  over 
6h  and  assume  that  the  PDF  for  0uf9i,  is  constant  for  angles  of  interest.  Then  Q'  becomes  independent  of  0t,  and  the  PDF  of  Q 
becomes 


MQ) 


n  j  d«jfajl«p  x(DSj)]dlQ  -  Q'K,  Dsj,  £)] 
Ylj  \  daj  fajLay,  MDsJ)\ 


(5) 


3.1.  The  General  Single  Screen  ( N  =  1 ) 
For  the  simple  case  of  a  single  scattering  screen,  equation  (5)  becomes 

,  (m  j  dafa{a-,-aD'^8[_Q  -  Q'(a,  Ds,  D)] 
2  j  dafa(a;—aD's) 

where 


(6) 


D's  =  DJt  1  -  DJD) . 


For  the  angle  of  arrival,  Q  =  0,Q'  =  a(DJD),  and 


f  m  =  i'dY  Jal(D/Ds)O-,-0(D  -  Dj] 
Je  )  \Dj  j  dafa(a;  —  aD's) 


(7) 


(8) 


For  the  TOA,  Q  =  t  and  Q'  =  Ds(  1  —  DJD)  \  a  |2/2 c,  so  only  the  magnitude,  |  a  \ ,  is  constrained.  Transformation  of  d(Q  —  Q') 
to  <5(  |  a  |  —  a, )  yields 


where  a.  is  a  unit  vector, 


f()  =  .  c\  djfja, a4> ;  —  a, a,,, D() 

JA>  \D'J  j  dafM-aDJ 


=  cos  4>x  +  sin  (j)y  , 


(9) 


(10) 


and 


(2  ct\112 


(11) 


The  flux  density  of  a  source  is  conserved  only  for  an  infinite  screen  with  homogeneous  statistics  because  only  in  that  case  is 
as  much  flux  scattered  toward  the  observer  as  is  scattered  away.  We  define  the  flux  dilution  factor  as  the  ratio  of  the 
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probability  that  rays  reach  the  observer  to  the  probability  for  a  uniform,  infinite  screen : 

n f  =  f  dafa(a;-aD's) , 


(12) 


equal  to  the  denominator  of  equation  (6).  For  a  uniform,  infinite  screen,  rjF  =  1.  In  general,  rjF  <  1. 

We  illustrate  these  expressions  by  considering  specific  cases. 

3.1.1.  Infinitely  Wide  Screen  with  Homogeneous  Statistics 

For  an  infinite  screen  with  homogeneous  statistics,  the  denominator  of  equation  (6)  is  unity.  Specializing  to  circularly 
symmetric /„,  we  find  the  normalized  one-dimensional  PDF  for  the  magnitude  9  =  \  0  \  : 


W)  =  2n9^jfJD0/Ds) . 

If /„  is  a  Gaussian  function  with  rms  angle  ca  in  each  coordinate  direction,  then 

m 


(13) 


xae29e-e2l2a% 

(14) 

=  (DJDK , 

(15) 

=  Tq  1e~tlX0U(t) , 

(16) 

=  D'sa2/c  . 

(17) 

The  one-dimensional  rms  angle 


=  ^=<|0|2>1/2 


(18) 


characterizes  the  observed  range  of  angles  and  in  this  case  is  identical  to  oe,  which  is  the  scaled  range  of  scattering  angles.  In 
general,  9rmi  a e. 

If  oa  oc  v  2  as  for  a  plasma,  then  the  AOA  PDF  has  a  scale  oe  oc  v  2  and  the  TOA  PDF  has  scale  r0  oc  v  4.  These  scaling 
laws  for  observable  quantities  rely  on  the  assumption  that  the  screen  is  infinitely  wide. 

3.1.2.  Circular  Screen  with  Finite  Radius 

Consider  a  circular  screen  with  radius  2fmax  centered  on  the  line  of  sight.  Now  the  PDF  for  9  is  truncated  for  9  >  0max  = 
Xm„/(D  —  Ds).  Again  adopting  circularly  symmetric,  Gaussian  statistics  for  a,  we  find 

m  =  K2(l  -  e-^T^e-^Uid^  -  9) ,  (19) 

C  =  \  ,  (20) 

where  the  unit  step  function  enforces  truncation  of  the  PDF  for  9  >  0max.  If  the  rms  scattering  angle  is  small,  <rn 
Xm,,,  D/DSD  —  Ds ),  the  scaling  of  the  observed  size  with  frequency  is  according  to  that  of  a a.  For  larger  scattering  angles,  the 
physical  size  of  the  screen  becomes  important.  To  see  this,  we  calculate  the  rms  angular  size,  which  is,  for  the  circularly 
symmetric  Gaussian  and  a  frequency  scaling  <t„  =  <ra0(v/v0)-2, 


0  _  ( Ds\{  v  \  Tl  -  (1  +  Qe  CT/2 

rms  \DJ\vJ  l  l-e-? 


(21) 


The  frequency  scaling  is  £  oc  v4.  At  large  v,  (  ^  1  and  9ri 
The  TOA  has  a  PDF  and  mean  value 


.  oc  v  2.  At  small  v,  £  ->  0,  and  9rms  becomes  frequency  independent. 


e~‘lz  o 

m  =  — - - F(tmax  -  t) ,  (22) 

r0(l  ~e 

<(>  =  „[l  -  {(y^)]  •  (23) 

fmax  =  •  (24) 

For  a  wide  screen,  £  ->  oo,  tmax  -*  oo,  and  <t)  =  r0,  as  before.  However,  a  narrow  screen  with  f  1  gives  <f>  ~  Cr0.  The  TOA 
scaling  is  <t>  oc  v-4  for  wide  screens  but  becomes  frequency  independent  for  very  narrow  screens. 
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The  observed  flux  of  a  source  viewed  through  the  truncated  screen  is  attenuated  by  the  scattering.  The  flux  dilution  factor  is 

nF  =  1  -  ■  (25) 

For  £  1,  qF  ~  £,  while  £  >  1  yields  rjF  =  1.  When  the  angular  diameter  and  pulse  broadening  of  a  source  are  observed  to 

have  anomalous  frequency  dependence,  i.e.,  when  scattering  is  dominated  by  a  truncated  screen,  the  flux  is  diminished. 
Inspection  of  equations  (21)  and  (25)  indicates  that  there  can  be  a  correlation  of  rms  angular  size  and  flux  density. 


3.2.  A  Screen  with  Arbitrary  Variations  of  Scattering  Strength 

Above,  we  considered  screens  with  extreme  variations  of  scattering  strength :  discontinuous  or  truncated  to  zero.  Here  we 
consider  other  cases  that  may  have  relevance  to  sources  that  are  viewed  through  scattering  regions  with  structure  on  scales 
~(D  —  Ds)9.  As  in  equation  (5),  we  specify  the  scattering  angle  a}  from  the  jth  screen  by  a  PDF  that  depends  on  location  along 
the  screen  x(Dsj) :  faj[_aj ;  x(Dsj)f  For  simplicity,  we  drop  the  j  subscripts  and  discuss  scattering  from  a  single  screen.  Also,  for 
ease  of  discussion,  we  consider  fa  to  be  a  two-dimensional  Gaussian  with  angular  variance  ajx)  that  varies  across  the  screen. 

In  some  circumstances,  the  angle-of-arrival  distribution  fe  is  determined  purely  by  the  shape  and  width  of fa,  while  in  others 
it  is  determined  by  the  variations  in  aa  across  the  screen.  If  oa  is  constant  across  the  screen,  f0  is  a  scaled  version  of /„  and  the 
frequency  scaling  of  the  observed  angular  size  is  identical  to  that  of  the  scattering  angles. 

However,  if  oa  varies  across  the  screen,  the  observed  angular  size  may  reflect  variations  of  aa  in  addition  to  or  rather  than 
the  width  of fa  itself.  Let  (a  be  the  characteristic  length  scale  on  which  oa  varies  across  the  screen.  We  compare  this  with  the 
observed  angular  diameter  projected  back  to  the  screen,  yielding  a  length  scale  { L  ~  (D  —  Ds)9rms.  We  consider  three  cases: 

1.  If  t?a  <  and  variations  in  a a  are  statistically  homogeneous,  the  variations  average  out.  The  scaling  with  frequency  of 
the  observed  AOA  width  is  identical  to  that  of  a,  which  is  due  to  the  microphysics.  Note  that  our  examples  of  truncated 
screens  in  previous  sections  satisfy  the  inequality  but  are  not  statistically  homogeneous. 

2.  If  f  a  P  ( L,  then  oa  «  constant  over  the  part  of  the  screen  sampled,  and  the  angle-of-arrival  PDF  is  determined  by  fa.  The 
frequency  scaling  is  again  determined  by  the  microphysics.  Depending  on  how  large  ta  is,  eventually  time  variations  in  the 
image  are  expected  on  a  timescale  ta  ~  fJVc ff,  where  Kff  is  an  effective  velocity  determined  by  the  velocities  of  the  source, 
medium,  and  observer  (e.g.,  Cordes  &  Rickett  1998). 

3.  If  ta  ~  the  AOA  PDF  is  determined  by  a  combination  of fa  and  the  spatial  variation  of  oa.  The  truncated  screen  of 
§  3.1.2  is  an  extreme  example  of  this  case.  The  frequency  scaling  is  likely  to  depart  significantly  from  that  of  the  microphysics. 


4.  EXAMPLES 

4.1.  Illustration  of  Anomalous  Frequency  Scaling 

Figure  1  shows  the  “ image  ”/«» the  pulse-broadening  function  f,  and  the  scaling  with  frequency  of  the  pulse-broadening 
time  for  two  cases:  (1)  an  infinite  screen  ( heavy  solid  lines )  and  (2)  a  truncated  circular  screen  centered  on  the  direct  line  of  sight 
( light  and  dashed  lines).  The  underlying  scattering  function  {fa )  is  a  circular  Gaussian  PDF,  and  the  rms  scattering  angle  scales 
as  v-2.  For  the  infinite  screen,  the  pulse  broadening  scales  as  v-4,  as  expected.  However,  the  truncated  screen  yields  truncated 
images  and  truncated  pulse-broadening  functions  if  the  rms  scattering  angle  is  large  enough  that  rays  from  the  screen’s  edges 
reach  the  observer.  Thus,  truncation  occurs  at  low  frequencies  and  not  at  high  frequencies.  For  the  example  given,  the  break 
frequency  is  approximately  0.5  GHz.  Actual  break  frequencies  will  depend  on  particular  sizes  and  scattering  strengths  of 
screens. 

As  a  second  illustration  of  anomalous  frequency  scaling,  we  calculate  the  pulse  broadening  for  a  pulse  scattered  by  a 
two-component  screen.  The  first,  circular  component  with  radius  X1  is  centered  on  the  line  of  sight  and  scatters  radiation 
much  more  strongly  than  the  remainder  of  the  screen,  which  is  of  infinite  extent.  The  distributions  of  AOA  and  TO  A  follow 
from  the  master  equations  (8)  and  (9).  The  rms  angular  size,  the  mean  TOA,  and  related  quantities  are 


„  /x  (D.'fl  +  MxXl  +  Ue-O-a  +  CJe-t'Ji* 

0rms  =  V2M,Z)JL  l+e-b-e-fc  J  * 

(26) 

/.N  T1  +  (Cl  +  uAi)e“?2  -  (1  +  Ci)e“Cll 

<(>  =  '{  i  +  *-<■  - J  • 

(27) 

U,2  =  , 

(28) 

Hfefc)'' 

(29) 

C2  =  -  Cl , 

(30) 

*2 


where  odl  2  is  the  rms  scattering  angle  produced  by  each  screen. 
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Fig.  1. — Scattering  from  an  infinite  screen  ( solid  lines )  and  a  truncated  circular  screen  of  radius  6  mas  (dashed  lines).  Top  left  panel :  Image  contours  for  an 
rms  angle  rre  =  10  mas.  For  the  infinite  screen,  the  contour  is  at  1  a  of  the  Gaussian  profile;  for  the  truncated  screen,  the  contour  represents  the  sharp  edge  of 
the  screen.  Top  right  panel :  Pulse-broadening  functions  on  a  loglinear  scale.  Middle  panel :  Mean  pulse-broadening  times  plotted  against  frequency  for  ae  =  1 
mas  at  1  GHz,  Ds  =  0.33  kpc,  and  D  =  1  kpc.  Bottom  panel :  Flux  dilution  factors  t]F. 


Figure  2  shows  0rms  plotted  against  frequency  for  different  ratios  adJad2  and  assuming  that  aiia  oc  v~2.  Note  that  we  vary 
odJod2  while  keeping  <rdl  constant.  The  figure  demonstrates  how  the  stronger  central  component  dominates  the  apparent 
source  size  at  high  frequencies  and  the  weaker,  distributed  component  dominates  at  low  frequencies.  At  intermediate 
frequencies,  there  is  a  plateau  where  the  angular  size  is  nearly  independent  of  frequency. 

Figure  3  shows  a  similar  plot,  now  for  the  mean  pulse-broadening  time,  <f>,  plotted  against  frequency  for  different  ratios 
t1(/t2  and  assuming  that  adia  °c  v-2.  The  roles  of  the  central  and  distributed  components  are  the  same  as  for  the  angular 
scattering  shown  in  Figure  2. 


4.2.  Scattering  from  Filaments 

Figure  4  shows  scattering  from  a  filament  located  along  the  direct  ray  path  for  three  values  of  rms  scattering  angle  in  the 
filament.  Small  rms  scattering  yields  a  circular  image  and  an  exponential  pulse-broadening  function.  For  sufficiently  large 
scattering  angles,  the  image  becomes  elongated  and  tends  toward  a  t~  1/2e-t  pulse-broadening  function 

Figure  5  shows  scattering  from  a  filament  at  different  locations  relative  to  the  direct  ray  path  but  for  identical  rms  scattering 
angles  in  the  filament.  When  the  filament  is  near  enough  to  the  image  center,  the  pulse-broadening  function  is  bimodal.  When 
the  filament  is  far,  the  pulse-broadening  function  is  dominated  by  the  much  weaker  scattering  from  outside  the  filament. 
Clearly,  a  superposition  of  filaments  near  the  direct  ray  would  produce  a  multiplicity  of  pulses. 

Figure  6  shows  scattering  from  an  ensemble  of  filaments  at  different  frequencies.  As  can  be  seen,  the  pulse-broadening 
function  shows  multiple  peaks  that  align  at  different  frequencies.  The  number  of  filaments  that  are  “  lit  up  ”  by  the  scattering 
decreases  in  going  to  higher  frequency.  The  scattering  screen  consists  of  a  very  strong,  extended  component  that  has 
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Fig.  2. — The  rms  angular  size  plotted  against  frequency  for  a  two-component  scattering  screen,  with  a  strongly  scattering  central  part  and  a  more  weakly 
scattering  extended  part.  The  curves  are  evaluations  of  equation  (26)  and  correspond  to  different  ratios,  adllod2  =  500,  50,  10,  5,  1  in  order  of  thickest  to 
thinnest  lines,  where  adl  is  the  rms  scattering  angle  from  the  central  part  and  ad2  is  the  value  for  the  extended  part.  In  producing  this  figure,  we  have  varied 
adiladi  while  keeping  adl  constant.  The  wavy  features  in  some  of  the  curves  are  real  and  reflect  the  different  frequency  dependences  of  terms  that  are 
subtracted  in  the  expression  for  rms  angular  size. 


embedded  “gaps”  where  the  scattering  is  weaker  but  still  strong  enough  to  scatter  radiation  toward  the  observer.  The 
extended  component  scatters  radiation  to  such  wide  angles  that  it  produces  negligible  contributions  to  the  pulse-broadening 
function  and  to  the  image.  Thus  the  strongest  contributions  to  measured  quantities  come  from  the  filamentary  gaps.  We  will 
explore  this  result  further  in  a  separate  paper  to  discuss  giant  pulses  from  the  Crab  pulsar. 

5.  APPLICATION  TO  GALACTIC  SCATTERING 

Galactic  sources  show  a  wide  range  of  scattering  levels,  indicative  of  the  concentration  of  intense  scattering  into  a 
Population  I-type  Galactic  distribution  (Cordes  et  al.  1985;  Taylor  &  Cordes  1993).  Here  we  discuss  particular  objects  whose 
scattering  may  be  interpreted  in  the  context  of  this  paper’s  formalism.  Our  discussion  is  brief.  We  defer  to  separate  articles  any 
detailed  treatment  on  particular  sources. 


5.1.  The  Crab  Pulsar 

The  Crab  pulsar  shows  enhanced  pulse  broadening  from  nebular  material  that  has  been  recognized  since  shortly  after  the 
pulsar’s  discovery  in  1968  (e.g.,  Vandenburg  1976  and  references  therein).  The  nebular  contribution  is  highly  episodic,  with 
dramatic  increases  of  the  pulse-broadening  time  by  a  factor  of  100  (Lyne  &  Thorne  1975;  Isaacman  &  Rankin  1977).  Recently, 
multiple  images  have  been  inferred  from  the  presence  of  echoes  of  the  pulse  shape  (Graham  Smith  &  Lyne  2000;  Backer  2000; 
Backer,  Wong,  &  Valanju  2000). 

Giant  pulses  from  the  Crab  pulsar  show  additional  evidence  for  nebular  contributions  to  the  scattering  that  are  probably 
from  discrete  filaments.  At  relatively  high  frequencies  (1.4-5  GHz),  giant  pulses  show  multiple  pulse  components  that  tend  to 
have  exponential  scattering  tails  with  time  constants  that  often  differ,  even  within  the  same  spin  period  of  the  pulsar  (Hankins 
&  Moffett  1998;  Sallmen  et  al.  1999;  Hankins  2000).  Also,  the  widths  of  the  pulse  components  appear  to  scale  less  strongly 
with  frequency  than  v-4.  These  characteristics  suggest  consistency  with  the  overall  picture  developed  in  this  paper.  A  detailed 
analysis  of  the  Crab  pulsar’s  pulses  is  deferred  to  another  paper. 
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Fig.  3. — Pulse-broadening  time  plotted  against  frequency  for  a  two-component  scattering  screen,  with  a  strongly  scattering  central  part  and  a  more 
weakly  scattering  extended  part.  The  curves  are  evaluations  of  equation  (27)  and  correspond  to  different  ratios,  iq/ij  =  500,  50, 10,  5, 1  in  order  of  thickest  to 
thinnest  lines,  where  t,  is  the  pulse  broadening  time  of  the  central  part  (if  it  were  infinite  in  extent)  and  t2  is  the  value  for  the  extended  part. 


5.2.  NGC6334B 

The  largest  angular  broadening  measured  is  for  the  extragalactic  source  NGC  6334B  viewed  through  the  H  ii  complex 
NGC  6334  (Trotter,  Moran,  &  Rodriguez  1998  and  references  therein),  9d  x  3"  at  1.4  GHz.  The  image’s  position  angle  rotates 
in  going  from  low  to  high  frequency,  and  the  image  axial  ratio  may  increase  from  1.2  at  1.4  and  5  GHz  to  1.5  at  15  GHz. 
Trotter  et  al.  (1998)  interpret  this  variation  as  signifying  an  outer  scale  for  the  wavenumber  spectrum  of  anisotropically 
scattering  density  irregularities,  *fout  <  1016  cm.  This  proposed  outer  scale  is  comparable  to  the  thickness  of  the  H  n  shell  in 
this  region  (Rodriguez,  Canto,  &  Moran  1988;  Kahn  &  Breitschwerdt  1990).  An  alternate  possibility  is  that  the  density 
irregularities  scatter  isotropically  but  the  anisotropic  images  reflect  density  irregularities  that  are  confined  to  the  H  n  shell. 

If  the  latter  were  the  case,  we  would  expect  that  the  axial  ratio  would  increase  as  a  function  of  decreasing  frequency, 
opposite  to  what  is  observed.  At  lower  frequencies,  the  size  of  the  minor  axis  of  the  scattering  diameter  would  be  constrained 
by  the  width  of  the  H  n  shell,  while  the  size  of  the  major  axis  would  be  essentially  unconstrained,  unless  the  scattering  is  so 
intense  that  the  size  of  the  major  axis  is  also  limited  by  the  scale  of  the  H  n  shell.  For  NGC  6334B,  the  axial  ratio  appears  to  be 
constant  with  frequency  or  increasing  with  increasing  frequency.  The  frequency  behavior  of  the  axial  ratio  indicates  that  any 
relevant  length  scales  in  the  H  n  shell  must  be  smaller  than  1014  cm  (the  smallest  length  scale  probed  by  the  highest  frequency 
observations)  or  larger  than  1017  cm  (the  largest  length  scale  probed  by  the  lowest  frequency  observations). 

5.3.  CygnusX-3 

The  compact  source  Cyg  X-3  is  heavily  scattered  (0'.'5  at  1  GHz)  and  has  an  anisotropic  image  (axial  ratio  ss 1.2  and 
increasing  with  increasing  frequency)  whose  position  angle  varies  with  frequency  (Wilkinson,  Narayan,  &  Spencer  1994; 
Molnar  et  al.  1995).  The  variation  of  position  angle  has  been  interpreted  by  Wilkinson  et  al.  (1994)  as  due  to  a  changing 
orientation  of  (anisotropic)  diffracting  irregularities  on  a  length  scale  of  order  9(D  —  Ds)  ~  0.01  pc.  In  their  picture,  the  image 
asymmetry  is  due  to  anisotropic  diffracting  irregularities,  and  the  orientation  change  with  frequency  is  attributed  to  the 
spatial  variation  of  those  irregularities. 

An  alternative  explanation  is  that  the  diffracting  irregularities  are  isotropic  and  the  image  anisotropy  reflects  spatial 
variations  of  the  strength  of  the  diffracting  irregularities  on  the  scale  9(D  —  Ds).  Molnar  et  al.  (1995)  have  proposed  that  the 
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Fig.  4. — Images  and  pulse-broadening  functions  for  a  vertical  filament  centered  on  the  direct  line  of  sight  but  with  different  rms  scattering  angles.  The 
filament  has  a  width  of  0.28  mas,  a  length  of  300  mas,  and  truncates  sharply.  Equations  (8)  and  (9)  were  evaluated  to  obtain  the  plots.  Left  panels :  Images 
shown  on  a  frame  of  size  2.5  x  2.5  mas.  The  +  symbol  designates  the  direct  line  of  sight  to  the  source.  The  rms  scattering  angle  is  constant  in  the  filament  and 
is  1  mas  v~2  =  1, 0.25,  and  0.11  mas  for  v  =  1, 2,  and  3  GHz  in  going  from  the  top  to  the  bottom  frame.  The  contour  levels  extend  from  the  maximum  down  to 
10-3  of  the  maximum  in  uniformly  spaced  intervals  of  0.5  decades.  Right  panels :  The  pulse-broadening  function  is  shown  on  a  loglinear  scale  with  1  decade 
on  the  vertical  axis.  The  horizontal  axis  extends  to  1.5  /is  in  the  top  frame  and  to  factors  of  16  and  81  times  smaller  in  the  middle  and  lowest  panels  to 
compensate  for  the  v-4  dependence  of  the  pulse-broadening  time  (if  there  were  a  continuous  screen).  The  vertical  dashed  line  indicates  zero  delay.  Top :  Large 
rms  scattering  angle  so  that  the  image  shape  is  dominated  by  the  edges  of  the  filament.  The  pulse-broadening  function  is  nearly  of  the  form  t~  1/2e~‘.  Middle : 
Smaller  rms  scattering  angle  so  that  the  image  is  less  elongated.  The  pulse-broadening  function  is  exponential  at  small  delays  but  shows  a  break  point  because 
there  is  an  absence  of  large  delays  due  to  the  truncation  of  the  filament.  Bottom-.  Rms  scattering  angle  small  enough  so  that  the  edges  of  the  filament  are  not 
seen.  The  pulse-broadening  function  is  exponential  in  form,  e-'. 


H  ii  region  DR  1 1  is  responsible  for  the  bulk  of  the  scattering  along  this  line  of  sight.  A  key  difficulty  with  this  explanation  is 
that  the  observed  axial  ratio  increases  with  increasing  frequency.  This  is  inconsistent  with  the  notion  of  scattering  from  a 
single  filament  but  may  be  consistent  with  scattering  from  a  group  of  filaments  that  are  individually  smaller  than  1016  cm. 

5.4.  Sgr  A*  and  Galactic  Center  OH/IR  Masers 

Galactic  center  sources  show  large  scattering  diameters  (~  1"  at  1  GHz)  and  significant  image  asymmetries  that  vary  across 
the  roughly  30'  size  of  the  scattering  region  (Lazio  &  Cordes  1998).  The  major  axis  of  the  image  of  Sgr  A*  shows  a  v~2 


1006 


CORDES  &  LAZIO 


Vol.  549 


Fig.  5. — Images  and  pulse-broadening  functions  for  a  filament  at  different  locations  from  the  direct  line  of  sight,  evaluated  using  equations  (8)  and  (9).  The 
filament  is  0.075  mas  wide  and  300  mas  long,  and  the  rms  scattering  angles  are  4  and  0.5  mas  inside  and  outside  the  filament,  respectively.  Left  panels :  Images 
of  size  2.5  x  2.5  mas  shown  with  contours  spaced  at  0.5  decade  intervals  down  to  10“ 3  from  the  peak.  The  +  symbol  designates  the  direct  line  of  sight  to  the 
source.  Right  panels'.  Pulse-broadening  functions  shown  on  loglinear  scales;  the  vertical  axis  covers  1  decade  and  the  horizontal  axis  covers  0.3  /is.  Top :  The 
image  is  a  circular  Gaussian  combined  with  a  subimage  from  the  filament.  The  filamentary  subimage  is  weaker  in  amplitude  because  flux  is  diluted  by  the 
heavier  scattering  from  the  filament.  The  pulse-broadening  function  is  nearly  exponential  in  form,  oce~‘.  Middle :  The  filament  is  close  enough  so  that  its 
subimage  now  contributes  significantly  to  the  pulse-broadening  function,  which  is  distorted  by  the  depletion  of  flux  at  some  delays  by  the  scattering  from  the 
filament.  Bottom :  The  circular  and  filamentary  subimages  are  sufficiently  merged  so  that  the  pulse-broadening  function  is  bimodal:  a  narrow  component 
from  the  portion  of  the  filament  nearest  the  direct  line  of  sight  and  a  second  peak  associated  with  the  nonfilamentary  scattering.  Relative  to  the  peak,  the 
pulse-broadening  function  is  now  slightly  broader  than  in  the  top  and  middle  cases. 


frequency  dependence  from  1.4  to  43  GHz  (Yusef-Zadeh  et  al.  1994;  Lo  et  al.  1998;  Bower  &  Backer  1998),  and  the  image  itself 
displays  no  change  in  its  major  axis,  axial  ratio,  or  orientation  on  timescales  as  long  as  a  decade  (Lo  et  al.  1998;  Marcaide  et 
al.  1999). 

Deviations  from  a  v  "  2  dependence  for  the  major  axis,  minor  axis,  or  both  are  claimed  at  a  variety  of  frequencies  from  as  low 
as  43  GHz  (Lo  et  al.  1998)  to  as  high  as  215  GHz  (Krichbaum  et  al.  1998).  These  deviations  from  the  v”  2  dependence  observed 
at  lower  frequencies  are  commonly  interpreted  as  an  effect  of  the  intrinsic  source  diameter  becoming  important  at  the  various 
frequencies. 

The  deviations  could  also  be  symptomatic  of  the  size  scale  of  the  scattering  medium.  If  this  is  the  case,  then  the  frequencies 
at  which  the  breaks  occur  correspond  to  spatial  scales  £  ~  5  AU  (for  a  break  occurring  at  43  GHz)  to  as  small  as  0.17  AU  (at 
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Fig.  6. — Images  and  pulse-broadening  functions  for  an  ensemble  of  filaments  superposed  with  a  much  stronger  scattering  screen.  The  filaments  therefore 
represent  gaps  where  the  scattering  is  weaker  than  in  the  screen  but  still  strong  enough  to  scatter  radiation  back  to  the  observer.  Physical  filament  widths  and 
lengths  are  0.03  and  1.5  mas,  respectively,  though  their  apparent  lengths  are  smaller  because  the  flux  scattered  to  the  observer  depends  on  distance  from  the 
direct  ray  path.  At  1  GHz,  the  screen  has  an  rms  scattering  angle  of  2  mas,  while  the  filaments  have  scattering  angles  assigned  randomly  and  uniformly 
between  0.1  and  2  mas.  Scattering  angles  scale  as  v-2.  From  top  to  bottom,  v  =  0.1,  0.2,  0.3,  and  0.5  GHz.  Left  panels:  Images  are  shown  in  2.5  x  2.5  mas 
frames  with  contours  spaced  at  0.5  decade  intervals  over  1  decade  from  the  image  peaks.  Contributions  from  the  strong  screen  are  too  small  to  be  reflected  in 
the  contours.  Right  panels:  Pulse-broadening  functions  shown  on  a  loglinear  scale  with  1  decade  on  the  vertical  axis.  The  horizontal  axis  covers  6.7  /is.  The 
observed  shapes  are  dominated  by  the  filamentary  gaps. 


215  GHz).  These  spatial  scales  are  comparable  to  the  outer  scale  inferred  by  Lazio  &  Cordes  (1998),  on  the  basis  of  a 
comparison  of  scattering  strength  and  thermal  free-free  emission.  However,  the  velocities  in  the  Galactic  center  (~  50  km  s“ 1) 
suggest  that  variations  in  the  orientation  of  the  image  of  Sgr  A*  should  be  seen  on  timescales  of  order  0.1  yr(//l  AU)(u/50  km 
s-1)-1.  As  mentioned  above,  these  are  not  seen.  This  implies  either  that  all  striations  in  the  medium  are  oriented  in  the  same 
direction  or  that  the  intrinsic  source  size  is  in  fact  important  at  high  frequencies. 

The  lack  of  variations  in  the  image  orientation  indicates  that  there  are  not  likely  to  be  striations  or  other  structure  in  the 
scattering  screen  on  scales  smaller  than  about  1  AU.  Variations  do  occur,  though,  on  much  larger  scales  of  order  15' 
corresponding  to  spatial  scales  of  roughly  25  pc,  the  separations  between  Sgr  A*,  the  various  OH  masers,  and  other  scattered 
sources  (Lazio  et  al.  1999).  Effects  from  the  spatially  limited  scattering  described  in  this  paper  are  unlikely  to  be  seen  with  the 
current  census  of  Galactic  center  sources.  The  narrowband  nature  of  OH  masers  means  that  the  frequency  dependence  of 
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their  scattering  diameters  cannot  be  measured,  and  Sgr  A*  itself  is  obscured  below  1  GHz  due  to  free-free  absorption  by 
Sgr  A.  Detection  of  additional  radio  transients  (e.g.,  Zhao  et  al.  1992)  or  radio  pulsars  (Cordes  &  Lazio  1997)  may  allow  such 
effects  to  be  detected  at  frequencies  below  1  GHz. 

5.5.  Extreme  Scattering  Events 

“Extreme  scattering  events”  (ESEs)  are  events  identified  in  the  light  curves  of  several  active  galactic  nuclei  (AGNs;  Fiedler 
et  al.  1987, 1994)  and  two  pulsars  (Cognard  et  al.  1993;  Maitia,  Lestrade,  &  Cognard  1999).  They  are  roughly  consistent  with 
refractive  defocusing  and  caustic  formation  from  discrete,  small-scale  plasma  structures  (Fiedler  et  al.  1987;  Romani  et  al. 
1987;  Clegg,  Fey,  &  Lazio  1998).  However,  alternative  explanations  invoke  the  outer  ionized  regions  of  predominantly 
neutral,  primordial  H  i  clouds  (Walker  &  Wardle  1998)  or  distributed  fluctuations  much  like  those  that  account  for  the 
diffractive  scintillations  of  pulsars  (Fiedler  et  al.  1994).  The  fundamental  difference  between  these  models  is  the  implied  gas 
pressure.  As  discussed  further  below,  a  discrete  structure  is  necessarily  overpressured  compared  to  the  general  ISM,  so  any 
such  structures  must  exist  either  transiently  or  in  regions  of  small-volume  filling  factor  that  can  support  such  pressure. 

If  ESEs  result  from  discrete  ionized  structures,  then  effects  described  in  this  paper  should  be  present  in  the  scattered  image 
and  pulse  shape  (for  pulsars).  This  notion  is  little  explored  because  few  ESEs  have  been  identified  and,  given  that  most  are  seen 
from  AGNs,  intrinsic  source  size  effects  can  also  diminish  the  appearance  of  diffraction  effects.  VLBI  observations  of  the 
source  1741  —  038  undergoing  an  ESE  have  shown  no  indication  of  a  truncated  image  (Lazio  et  al.  2000).  However,  those 
observations  were  at  relatively  high  frequencies  (>1.7  GHz)  and  had  limited  dynamic  range.  Pulse-timing  observations  of 
PSR  B1937  +  21  (Cognard  et  al.  1993)  show  no  change  in  the  pulse  shape,  though,  again,  these  observations  were  obtained  at 
1.4  GHz.  Future  observations  of  a  source  undergoing  an  ESE  at  lower  frequencies  (e.g.,  0.33,  0.41,  or  0.61  GHz)  would  place 
much  more  stringent  constraints  on  the  notion  that  ESEs  arise  from  discrete  ionized  structures. 

6.  FUTURE  OBSERVATIONS 

6.1.  Application  to  Extragalactic  Scattering 

Scattering  from  extragalactic  plasma  can  arise  from  the  distributed  intergalactic  medium  (IGM),  most  of  which  is  expected 
to  be  ionized,  from  intervening  Lya  clouds,  and  from  intervening  galaxies.  Of  greatest  relevance  to  this  paper  are  the  last  two 
cases,  and,  of  these,  intervening  galaxies  are  likely  to  be  the  more  important  because  of  their  greater  column  densities.  A 
face-on  galaxy  like  the  Milky  Way  will  scatter  radiation  from  a  background  source  into  an  apparent  size  of  at  least 
dd  ~  l(v/0.33  GHz)- 2  2  mas  (Cordes  &  Lazio  1991;  Taylor  &  Cordes  1993).  Scattering  by  H  ii  regions  yields  even  larger 
angles,  so  some  background  sources,  albeit  at  low-probability  alignments,  will  display  images  that  reflect  the  sizes  of  H  ii 
regions  and,  in  some  instances,  spiral  arms  that  contain  them. 

Scattering  from  an  edge-on  galaxy  will  be  about  103-104  times  larger,  or  1"— 10"  at  0.3  GHz.3  The  lateral  scale  is 
9dD  ~  (15-150)D3OOO  kpc  for  D  =  3000D3OOO  Mpc.  Thus,  near  edge-on  galaxies  will  produce  scattered  images  that,  in  part, 
display  the  shapes  of  the  galaxies.  At  even  lower  frequencies,  scattering  diameters  from  Lya  clouds  and  galaxies  with 
significantly  smaller  scattering  strength  will  produce  similar  effects.  Imaging  radio  observations  at  ~  0. 1  GHz  will  thus  probe 
intergalactic  structures. 

We  defer  to  another  paper  a  thorough  discussion  of  intergalactic  scattering,  taking  into  account  cosmological  expansion 
and  evolution.  Scattering  may  be  able  to  probe  the  intergalactic  medium  at  redshifts  near  the  reionization  epoch. 

6.2.  Low-Frequency  Galactic  Observations 

The  (nominally)  strong  frequency  dependence  of  interstellar  scattering  observables  suggests  that  the  anomalous  scattering 
described  here  will  most  likely  occur  at  low  frequencies.  High-resolution,  low-frequency  instruments  such  as  the  Giant 
Meterwave  Radio  Telescope  (Ananthakrishnan  1995)  and  the  proposed  Low-Frequency  Array  (Kassim  et  al.  2000)  and  the 
low-frequency  Square  Kilometer  Array  (Butcher  2000)  have  or  will  have  subarcminute  resolution  at  frequencies  below  150 
MHz.  Consequently,  they  may  detect  anomalous  scattering  along  lines  of  sight  less  heavily  scattered  than  those  described  in 
§  5.  Here  we  consider  relevant  lines  of  sight  and  frequencies  for  which  anomalous  scattering  is  a  possibility. 

The  relevant  length  scale  in  regions  of  less  intense  scattering  may  be  the  outer  scale  of  the  density  fluctuation  spectrum  *f0. 
(This  may  also  be  the  relevant  length  scale  in  intense  regions,  though  its  value  could  be  quite  different  and  potentially  much 
smaller.)  Near  the  Sun  (within  « 1  kpc),  l0  ~  1  pc  (Armstrong,  Rickett,  &  Spangler  1995).  It  is  unlikely  that  scattering 
diameters  will  probe  this  spatial  scale  (i.e.,  0dD  ~  Z0)  unless  v  <  10  MHz.  As  the  ionosphere  becomes  increasingly  opaque  at 
frequencies  v  <  10  MHz,  ground-based  interferometric  arrays  will  likely  not  be  affected  by  anomalous  scattering  in  the  solar 
neighborhood. 

Toward  the  inner  Galaxy,  Galactic  latitudes  | /|  <  50°,  stronger  scattering  than  that  in  the  solar  neighborhood  (but  still 
weaker  than  the  intense  scattering  described  in  §  5)  will  obtain.  In  this  case,  observations  at  meter  wavelengths  may  display 
anomalous  scattering. 

A  competing  effect  for  detecting  anomalous  scattering  is  free-free  absorption.  The  density  fluctuations  responsible  for 
interstellar  scattering  also  contribute  to  free-free  absorption.  Sources  seen  along  heavily  scattered  lines  of  sight  at  shorter 
wavelengths  may  be  free-free  absorbed  at  longer  wavelengths.  For  instance,  free-free  absorption  renders  the  Galactic  center 
increasingly  opaque  for  frequencies  v  <  1  GHz  (Anantharamaiah  et  al.  1991). 


3  This  large  increase  in  the  scattering  diameter  occurs  because  of  the  presence  of  enhanced  scattering  regions,  such  as  those  described  in  §  5.  In  the  Milky 
Way,  the  approximate  mean  free  path  between  enhanced  scattering  regions  is  8  kpc  (Taylor  &  Cordes  1993).  A  line  of  sight  through  the  disk  of  a  galaxy  like 
the  Milky  Way  is  quite  likely  to  encounter  one  or  more  enhanced  scattering  regions,  leading  to  the  large  increase  cited. 
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7.  IMPLICATIONS  FOR  THE  INTERSTELLAR  MEDIUM 

As  alluded  to  before,  the  existence  of  compact,  turbulence-containing  ionized  structures  is  directly  related  to  their  longevity 
and  rarity,  or  filling  factor,  in  the  Galaxy.  Except  for  chance  fluctuations  from  distributed  turbulence,  many  observed 
phenomena  suggest  the  existence  of  compact  structures  with  densities  that  imply  they  are  overpressured  compared  to  most  of 
the  ISM.  This  is  not  overly  surprising  because  the  ISM  is  highly  dynamic  and  is  in  pressure  equilibrium  only  in  some  average 
sense.  It  is  not  known  which  kinds  of  locales  (H  n  regions,  supernova  shocks,  etc.)  provide  the  largest  scattering  strengths.  We 
suggest,  simply,  that  the  observable  effects  described  in  this  paper  might  be  used  to  better  probe  the  physical  sizes  of  regions 
that  produce  the  largest  levels  of  scattering. 

Striations  in  interstellar  gas  densities  on  subparsec  scales  are  most  likely  associated  with  magnetic  fields.  On  diffraction 
scales  <  1011  cm,  a  compelling  idea  is  that  turbulence  is  essentially  two  dimensional  and  irregularities  are  elongated  along  the 
field  lines  (Higdon  1984, 1986;  Goldreich  &  Sridhar  1995).  Larger  scale  filaments,  such  as  those  seen  near  the  Galactic  center 
perpendicular  to  the  plane  of  the  Galaxy  (Yusef-Zadeh  &  Morris  1987),  are  also  thought  to  be  along  the  local  field  direction.  If 
anisotropically  diffracting  irregularities  are  contained  in  screens  that  are  themselves  elongated  in  the  same  direction,  it  may  be 
difficult  to  separate  (and  thus  identify)  the  two  possible  contributions  to  image  elongation. 

The  frequency  dependence  of  anomalous  scattering  may  offer  a  means  for  identifying  the  cause  of  image  elongation  for  a 
particular  source  if  it  is  heavily  scattered  by  a  single  (or  few)  filament(s)  (see  Fig.  2  and  §  5).  If  the  image  elongation  arises  from 
anisotropic  scattering  by  small-scale  density  irregularities,  increased  spatial  resolution  of  the  scattering  material  (e.g.,  by 
making  observations  at  higher  frequencies)  may  yield  increased  axial  ratios  as  less  spatial  averaging  is  done  over  the 
small-scale  irregularities  (e.g.,  as  argued  by  Wilkinson  et  al.  1994  for  the  image  of  Cyg  X-3).  This  change  in  axial  ratio  is 
expected  only  if  there  are  spatial  variations  of  the  orientation  of  scattering  turbules.  Conversely,  if  image  elongation  is 
produced  by  the  boundaries  of  filaments,  a  different  frequency  dependence  may  be  seen.  A  more  complicated  frequency 
dependence  may  result  if  the  scattering  results  from  a  number  of  smaller  filaments  (e.g.,  Fig.  6). 

The  wavenumber  spectrum  of  electron  density  irregularities  is  often  constrained  by  the  scaling  law  of  angular  size  and  pulse 
broadening  (and  its  reciprocal,  the  scintillation  bandwidth)  with  frequency.  For  moderate  scattering,  where  the  dominant 
length  scales  are  between  the  inner  and  outer  scales,  the  pulse-broadening  time  scales  as  v~x  with  x  =  —  2 /?/(/?  —  2),  where  fi  is 
the  exponent  of  the  three-dimensional,  isotropic  wavenumber  spectrum  (Cordes  et  al.  1985;  Rickett  1990).  For  a  Kolmogorov 
spectrum,  [S  =  11/3  and  *  =  4.4.  If  irregularities  are  isotropic,  so  that  [I  =  11/3,  but  the  medium  is  confined  in  the  transverse 
direction,  the  actual  value  of  x  is  lessened.  The  value  of  fi  inferred  would  be  greater  than  11/3  in  this  instance.  A  similar  trend 
occurs  when  angular  broadening  is  used  to  infer  (L  It  is  not  clear  which,  if  any,  of  the  published  constraints  on  /I  are  affected  by 
the  influence  of  scattering-region  confinement.  A  detailed  study  of  the  wavenumber  spectrum  is  deferred  to  another  paper. 

8.  SUMMARY 

In  this  paper  we  have  shown  that  radio  scattering  observables  such  as  image  shapes  and  pulse-broadening  functions  can  be 
strongly  influenced  by  structure  in  the  scattering  medium  on  length  scales  substantially  larger  than  those  that  cause  the 
scattering.  As  such,  careful  multifrequency  observations  can  be  used  to  constrain  properties  of  the  interstellar  medium  on 
scales  that  are  typically  ~  1-10  AU.  Intergalactic  scattering  has  not  been  identified  but  is  certainly  expected  from  intervening 
spiral  galaxies,  probably  expected  for  some  Lya  clouds,  and  may  occur  from  distributed  ionized  gas.  For  intergalactic 
scattering,  relevant  length  scales  can  be  comparable  to  the  sizes  of  galaxies.  A  low-frequency  VLBI  survey  of  extragalactic 
sources  may  thus  probe  the  level  of  scattering  in  other  galaxies  and  in  the  general  intergalactic  medium.  It  is  also  expected 
that  scattering  of  radiation  from  gamma-ray  burst  afterglows  will  be  influenced  in  some  cases  by  intervening  ionized  gas  in  the 
IGM  as  well  as  in  the  Milky  Way’s  ISM.  These  issues  will  be  explored  in  separate  articles. 
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